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Abstract 

The question of optimal portfolio is addressed. The conventional Markowitz portfolio optimi- 
sation is discussed and the shortcomings due to non-Gaussian security returns are outlined. A 
method is proposed to minimise the likelihood of extreme non-Gaussian drawdowns of the port- 
folio value. The theory is called Leptokurtic, because it minimises the effects from "fat tails" of 
returns. The leptokurtic portfolio theory provides an optimal portfolio for investors, who define 
their risk-aversion as unwillingness to experience sharp drawdowns in asset prices. Two types of 
risks in asset returns are defined: a fluctuation risk, that has Gaussian distribution, and a draw- 
down risk, that deals with distribution tails. These risks are quantitatively measured by defining 
the "noise kernel" - an ellipsoidal cloud of points in the space of asset returns. The size of the 
ellipse is controlled with the threshold parameter: the larger the threshold parameter, the larger 
return are accepted for investors as normal fluctuations. The return vectors falling into the kernel 
are used for calculation of fluctuation risk. Analogously, the data points falling outside the kernel 
are used for the calculation of drawdown risks. As a result the portfolio optimisation problem 
becomes three-dimensional: in addition to the return, there are two types of risks involved. Opti- 
mal portfolio for drawdown-averse investors is the portfolio minimising variance outside the noise 
kernel. The theory has been tested with MSCI North America, Europe and Pacific total return 
stock indices. 

PACS numbers: 89.65.Gh, 89.75.Da, 05.40.Fb, 05.45.Tp 
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INTRODUCTION 

The bridges between statistical physics and financial economics have been recently crossed 
by many of authors, cf. jll 0,0,0). The topics of research include the descriptive statistics 

s and many others. Recently the prob- 



!. The portfolio optimisation has 



of the price movements, market microstructure mode 
lem of basket of assets has attracted attention (cf. 
been introduced in 1950s by Harry Markowitz [6] . The simplicity of portfolio optimisation 
problem has made it well-accepted in financial community (cf. |7[ and references therein). 
The conventional or Markowitz portfolio theory (MPT) assumes the Gaussian probability 
distribution function for security returns - this is widely questioned in Econophysics liter- 
ature. The recent reports show that even the Levy stable distributions are not describing 
the stochastic process of price changes 0]. Furthermore, the temporal organisation of the 
price increments is also complicated, multifractal, cf. Apparently, the optimal portfolio 
question should be reconsidered bearing in mind these very important findings. However, 
these revised methods have to be simple and robust. Indeed, the statistical data available for 
the analysis are typically insufficient for a statistically meaningful application of advanced 
and complicated techniques, 

In this paper the theory is provided that leads to three-dimensional portfolio optimisation. 
In order to keep the approach as simple as possible, we ignore the multifractal aspects of 
the price movements. The risk of the portfolio is split to two: Gaussian risks and "fat 
tail" risks. Then, the portfolio choice is provided for the investors who want to minimise 
the absolute drawdowns. In this paper the Leptokurtic Portfolio Theory is introduced and 
motivated. Further, the application is provided. Finally, the approach is tested based on 
various international stock indices. 



Definitions 

Let rj(t, At) be the return of the security % in the portfolio at time t for recent period of 
At: 

r,(f, At) ee ln Pi (t) - ln Pi (t - At) ~ ^'^'^ W 

where Pi(t) denotes the price of security % at time t. As noted earlier, the statistical properties 
of quantity r(t) are well-elaborated in the literature. MPT assumes the Gaussian distribution 
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- an obvious disagreement with recent findings [8] . The portfolio return - simple sum of the 
returns of portfolio constituents and portfolio risk - standard deviation a p - are defined as 
follows: 

n 

r p (t,At) = J2 w Mt,At) (2) 

n n 

a l = Yl Yl W i W 3 a i,j ( 3 ) 
i=l j=l 

where cr^j is the covariation of returns of securities % and j; cr^ = of. Under the assumption 
of Gaussian returns, the measured historic values of returns and covariations can be used as 
a proxy to the future. 

The expected (i.e. ex ante) return in time-window r is defined as average of previous 
realisations: 

tft = E{r t ) = (r«> T (4) 
The covariance of returns for security pair is found as follows: 

°ij = ( r irj) T - (n) T (r j ) r = (r i r j ) r - ruVj ( 5 ) 

If % — j then the covariance relaxes to simple variance: = af = (rf) T — r/f. 

It is convenient to represent the portfolio in the matrix form. Denoting the return and 
weight vectors respectively with r\ and w, Eq. (J2J) is rewritten as 

r p = r] T w, (6) 

where 7] T is transpose of rj. Further, let C be the matrix of covariation coefficients with 
elements of ai t j. The Eq. (J3J) becomes: 

al = w T Cw (7) 



LEPTOKURTIC PORTFOLIO THEORY 



What is the quantitative measure of a portfolio risk in the case of strongly non-Gaussian 
price fluctuations with (possibly) infinite variance and Levy-like spectrum of price jumps? 
Note that stable Levy distribution of price jumps assumes diverging variance; the power 
law reported by X. Gabaix et al [8| assumes diverging kurtosis (fourth moment). In the 
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case of diverging fourth moment, the distributions are called leptokurtic. The tails of the 
leptokurtic distributions are "fatter" than predicted by Gaussian distribution. In fact, all 
power-law distributions lead to leptokurtic distributions and that must be accounted by 
determination of risk. The obvious question arises: what is the risk of portfolios that obey 
leptokurtic distribution tails? Here, a model is proposed to offer the portfolio choice under 
such assumptions. First, the two types of fluctuations are separated: the fluctuation or 
Gaussian risk (which can be called a "good" risk, because the distribution tails approach 
quickly zero) and the drawdown, or power-law risk ("bad" risk, because the actual portfolio 
loss can be much larger than that of predicted by Gaussian approach). These two types of 
risk, together with the expected return, create a three-dimensional space; second subsection 
is devoted to the discussion of the portfolio optimisation in that space. Finally, a simple 
empirical illustration is provided using international stock indices. 

Separation of risks 

The power-law distribution of security returns leads to a large amount of "small" price 
movements, and few "large" movements. It is convenient to separate the "noise kernel" of 
daily returns — these are the returns which are smaller than few standard deviations. First, 
let us assume that the volatilities are not diverging. Then, the "noise kernel" is defined as a 
set of points Q in the n- dimensional space of returns (where i-th axes measures the return 
Ti of the i-th asset) falling into such an ellipsoid where the probability density function of 
return vectors is above a (small) threshold. This threshold serves as a model parameter (see 
below). Using the eigenvectors of the covariance matrix as the orthonormal basis p«, the 
ellipsoid is defined as 

n 

QeXXAT^n*, (8) 
i=i 

where \ is the i-th eigenvalue of the covariance matrix. The threshold parameter ~ 1 
should be chosen in such a way that the probability for a return vector falling outside the 
ellipsoid, is of the order of few percents; it regulates the ellipsoid " size" in units of standard 
deviations. In the case of two assets, the ellipsoid becomes ellipse, and can be expressed as 

Q E r]a) + f)a] - 2r» M < [2(o?oJ - o*,)] 9 (9) 



4 



Hereinafter, f denotes the zero-shifted return, i.e. fj = — rji. The shift is necessary to keep 
the centre of the ellipse in the 0-point of the coordinates. 

For diverging volatilities, the above outlined approach can be still applied; however, the 
net volatility is then a non-stationary quantity and depends on the particular observed set 
of extreme price movements. The effect of these extreme movements can be removed by 
an iterative approach: instead of using the net covariance matrix for the calculation of the 
noise kernel, the noise kernel's own covariance matrix must be used. 

The investors are looking for capital appreciation, in order to achieve the investment 
performance targets. What the investors do not like, is fluctuation in asset prices. Indeed, 
everybody would be happy to see a linear appreciation of capital with a suitable growth rate. 
However, investors would accept the "normal" fluctuations with a low standard deviations. 
What the investors really want to avoid, is a sharp drawdown in asset prices. Thus, the risk 
components should be separated to Gaussian ("good") fluctuation risk, and drawdown risk. 
This is done as follows: 

Definition 1. Fluctuation risk: The risk is measured as the average deviation from 
average return, i.e. risk is measured with standard deviation. 

Definition 2. Drawdown risk: The risk is measured as the likelihood of sharp changes 
in asset prices. Risk is the degree of willingness to experience sharp moves in asset prices. 

According to MPT, the drawdown risk in not existing. This is actually not true, if one 
thinks back to the history: stock market crashes have occurred quite frequently: 1929 and 
1987, 2000 in US, 1997 South-East Asia, 1998 Russia are just some well-known examples. In 
order to separate the risks, the measure of correlation should be revised. Here, the concept 
of exceedance correlation (cf. jsl licl 111) is extended. The idea of separation of risks was 



also suggested by Chow et al (cf. [12J). 

Our simplified approach is as follows. For the calculation of the fluctuation risk of the 
portfolio (Def. only these data points which are in the noise kernel are accounted for. 
On the other hand, the points outside noise kernel are used for the calculation of drawdown 
risk of the portfolio, and of the related covariance matrix of returns (Def. EI). Note that 
the larger the threshold parameter 9, the larger fluctuations are considered as "normal" or 
"acceptable". To conclude, equations 
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%|fluc(0) = ^ f j)re Q ~ ( f i)reQ ( f i>reQ ( 10 ) 

and 

cr ij|draw(^) = ^'i)r3Q ~ (^)rsQ (^i)r3Q (H) 
are used to calculate the covariation coefficients for fluctuation risk and drawdown risk 
respectively. Here, (• • -) re n (and (• • -) r3 g) denote averaging over all the return vectors r 
belonging (and not belonging) to the kernel Q. The covariation matrix C used in MPT 
[cf. Eq (JJ|)] is thus split into C g Q and C 9 q, corresponding to the fluctuation and drawdown 
regimes, respectively. 

The separation of noise kernel is shown schematically in Fig. where the data points 
and their linear regression lines are plotted. In the left panel of Fig. ^ two normalised 

FIG. 1: Exclusion of the noise kernel from covariation analysis, a) Two normalised, but correlated 
(with p ~ 70%) Gaussian variables; b) daily returns of Standard & Poors's 500 and DAX equity 
indices. 




♦ KERNEL ■ Non-KERNEL 
Linear (KERNEL) Linear (Non-KERNEL) 



(with zero mean and unit standard deviation), but correlated (p ~ 70%) random Gaussian 
variables are plotted. The noise-kernel corresponds to the threshold parameter 9=1. We 
are dealing with a pure Gaussian distribution; therefore, the correlation does not change 
with the exclusion of the noise kernel. A different behaviour is observed in the right panel 
of Fig. ^ There are daily returns plotted for US Standard & Poors' 500 (spx) and German 
Deutsche Aktieindex (dax) equity indices. The period ranged from September 1959 to March 
2004 (i.e. total more than If 000 data points) and the kernel was defined with parameter 
9 = 3. As seen, the noise kernel and the rest of the data have quite different correlations. 



6 



Application of LPT 

Leptokurtic Portfolio Theory (LPT) is a simple generalisation of MPT. Equations (JHJ) 
and (Jll|) yield non-trivial results only when large asset price jumps (drawdowns) exist. The 
boundary between the ordinary fluctuation risk and extraordinary drawdown risk is defined 
by the parameter 9. Here, we have presented market data analysis with 6 = 3; the value 
6 = 2 yields similar results. Note that in the case of Gaussian distribution, the events 
fall into three standard deviations with probability equal to 99.73%, which seems to be a 
reasonable crossover point. 

The application of LPT introduces additional dimension to the Gaussian risk return space: 
non-Gaussian risks. So, the portfolios are evaluated in three-dimensional space containing 
two types of risks, and the return. This allows us to separate the fluctuation and drawdown 
risks and the portfolio choice gets more complicated. Investors must choose an optimal 
balance between the Gaussian risk and return, but they have also to evaluate the potential 
of the "out-of-statistics" drawdowns. To elaborate this concept, consider the example of 
Fig. db). Suppose the investor is willing to invest into a combination of the indices of SPX 
and DAX. If the investor is accepting the fluctuations within three standard deviations, but 
he is unwilling to accept the large drawdowns, then he should invest into a portfolio with 
minimised drawdown risk. On the other hand, if the low values of the average short-term 
Gaussian fluctuation is more highly prioritised than the desire to minimise the likelihood of 
drawdowns, minimised kernel-risk portfolio satisfies the investment goal. 

Equation provides the definition of returns in MPT. In the spirit of LPT, one could 
wish to distinguish between "ordinary" (Gaussian) and "extraordinary" returns. However, 
in the case of predicted returns, such splitting seems unjustified: first of all, because the 
use of historic returns as a proxy to the future returns is far from being justified. Indeed, 
due to the non-stable nature and non-stationarity of returns, there will be large differences 
between short-term realisations. This will inevitably lead to the fact that historic returns 
do not provide the valid forecast to the future returns - a circumstance that is easy to check 
in any financial time series. In fact, the expected return of the security depends also on 
the choice of investment horizon - it varies from investor to investor. To conclude, we leave 
the expected return un-splitted and prefer more fundamental (i.e. economic result based) 
approach in estimating the future returns. 
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Empirical evidence: the validity of LPT 

In previous, we have defined the idea of separation of risks and their measures. Here, the 
empirical evidence is provided. We construct the portfolios by using monthly data of the 
following Morgan Stanley Capital International (MSCI) total return indices: MSCI Europe 
(hereinafter denoted as E), MSCI North America (NA) and MSCI Pacific (P) in the period 
from December 1969 to February 2005 (total more than 420 months) . In this paper, we only 
look at the examples of pairs of securities. The portfolio problem with N securities will be 
addressed in future. The method used is the following: 

1. we take the pair of securities and find the noise kernels with 8 = 1,2,3 according to 
Eq. ©, 

2. we determine the kernel and non-kernel covariations according to Eqs. (|10j) and ([lljh 

3. we determine the minimum-risk portfolios according to fluctuation, drawdown and 
Markowitz definition (the latter is simply aggregated standard deviation of kernel and 
non-kernel) with parameter 8 = 1,2,3, 

4. we back-test the minimum-risk portfolios with the same data-set by investing hypo- 
thetical unit currency in December 1969 

This process is illustrated with an example of MSCI North America and MSCI Pacific 
time series. In Table [Q the correlation matrices are presented for MPT, Fluctuation, and 
Drawdown methods using the value 8 = 3. In Fig. the portfolio choice is shown using the 
same data as in Table HI The risk-return space in Fig. |2]is constructed in a spirit of MPT: 
in the abscissa, the risk in annualised standard deviations is plotted. Note that for different 
portfolio sets, the different definitions of covariations and standard deviations are used. The 
portfolios with minimum variance are marked with circles. 

TABLE I: Covariation matrices of MSCI North America and MSCI Pacific covariations 





MPT 


Fluctuation 


Drawdown 


&NA,P 


NA P 


NA P 


NA P 


NA 


0.1959% 0.1037% 


0.1365% 0.0797% 


0.8648% 0.3427% 


P 


0.1037% 0.3563% 


0.0797% 0.2647% 


0.3427% 1.4603% 



s 



FIG. 2: Determination of minimum risk portfolios of MS CI North America and MS CI Pacific time 
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Finally, a simple back-test is performed with portfolios found as follows: a hypothetical 
unit of currency is invested into optimal portfolios found in previous step as of December 
1969. The portfolio is re-balanced with initial allocations found on monthly basis. LPT 
is valid when the portfolios which are optimised for drawdowns provide smaller extreme 
drawdowns than portfolios which are optimised for variance (i.e. MPT portfolios). 

This test was carried out for all three above-mentioned security pairs for 9 e [1,5], total 
15 times. In Table ITTl the back-test is presented with previous example of NA and P. 



TABLE II: The return statistics of MSCI North America / Pacific portfolios that are minimized 
under different definitions of risk 



Risk type 


MPT risk 


Fluctuation risk 


Drawdown risk 


Min 


-23.3% 


-23.4% 


-23.13% 


Max 


14.39% 


14.1% 


14.92% 


Mean 


0.9% 


0.9% 


0.91% 


Standard deviation 


4.2% 


4.21% 


4.18% 


Mean (annualised) 


10.83% 


10.79% 


10.89% 


Std. dev. (annualised) 


14.54% 


14.59% 


14.48% 



Our results confirmed the theory: in all 15 cases the portfolio found by using LPT 
provided smaller absolute drawdown. Interestingly, all of the LPT portfolios had also higher 
return compared to those optimised for MPT. 
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DISCUSSION 



A simple and robust method of non-Gaussian portfolio optimisation has been devised. 
The concept of Gaussian noise kernel for the vector of daily asset return vectors led us to the 
separation of two types of risks. The portfolios that were designed to minimise drawdown 
risk provided higher average return with lower maximum drawdowns, than the minimum-risk 
portfolio according to the MPT. Although we have found very good results with the given 
data-sets, there is a clear need for further extensive analysis with real data and portfolios 
consisting of more than two assets. 

The support of Estonian SF grant No. 5036 is acknowledged. Robert Kitt also wishes to 
thank Hansa Investment Funds Ltd. 
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